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ABSTRACT 

Although the Eddington limit has originally been derived for stars, recently its rele- 
vance for the evolution of accretion discs has been realized. We discuss the question 
whether the classical Eddington limit - which has been applied globally for almost 
all calculations on accretion discs - is a good approximation if applied locally in the 
disc. For this purpose, a critical accretion rate corresponding to this type of modified 
classical Eddington limit is calculated from thin a-disc models and slim disc models. 
We account for the non-spherical symmetry of the disc models by computing the local 
upper limits on the accretion rate from vertical and radial force equilibria separately. 

It is shown that the results can differ considerably from the classical (global) 
value: The vertical radiation force limits the maximum accretion rate in the inner disc 
region to much less than the classical Eddington value in thin a-discs, while it allows 
for significantly higher accretion rates in slim discs. We discuss the implications of 
these results for the evolution of accretion discs and their central objects. 
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1 INTRODUCTION 

Similar to the stellar case, disc accretion may be limited by 
radiation pressure, counteracting gravity and viscous dis- 
sipation. The major difference to the stellar case is the 
different geometry which potentially complicates the situ- 
ation considerably. While in the stellar case, we are dealing 
with a - for all practical purposes - spherically symmetric, 
i.e., 1-dimensional situation, discs require an - at least - 
2-dimensional treatment. As yet, however, usually a quasi- 
stellar Eddington limit is being applied for disc models: (1) 
spherical symmetry of the system; (2) isotropic radiation; (3) 
homogeneous degree of ionisation; and (4) no time depen- 
dence (stationarity) . All these approximations do not apply 
for the disc case, though it is not clear a priori to what 
degree a proper treatment will alter the resulting numbers. 
In addition, independent of whether we are dealing with 
the stellar or the disc case, classically, the Eddington limit 
uses several approximations which may or may not be jus- 
tified in all cases: (i) Thomson scattering as the sole source 
of opacity; (ii) negligible gas pressure, in comparison with 
radiative pressure; and (iii) no relativistic effects. It is the 
goal of the present investigation to find out the importance 
of relinquishing some or all of these approximations. 

Over the at least two decades, a number of ef- 
forts have been carried out to investigate the ap- 
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plicab ility of an Eddington-type limit t o the disc 
case: [ J aroszyriski, Abramowicz fc Pacz vnskil il980tl and 
lAbramowicz. Calvani fc Nobilil il980fl . for instance, have 
been working on supercritical accretion discs. As a main re- 
sult of their work, models for so-called slim and thick accre- 
tion discs have been developed. In addition to standard ra- 
diative cooling, in these discs also cooling by advective flows 
is taken into account. This permits the critical accretion rate 
to increase without rising the radiative flux accordingly. The 
photo n trapping mechanism feegelmanlll978t lOhsuea et alJ 
2002) amplifies this effect even more. 

ICollin fc Kawaeuchil l|2004h inspected Narrow Line 
Seyfert Galaxies 1 and concluded that these accrete at super- 
Eddington rates, while their luminosity stays of the order 
of the classical Eddington limit. The authors proposed an 
additional non-viscous energy release in the gravitation- 
ally unstable region of the disc, which emits a frac tion of 
the optical luminosity. Recently. iTurner et al.l i2005t) inves- 
tigated the effect of photon bubble instabilities by two- 
and three-dimensional numerical radiation MHD calcula- 
tions, finding that photon bubbles may be important in cool- 
ing radiation-dominated accretion discs a nd therefore con- 
fir ming previous work by lGjnrmiie] il998T) , iBegelmanl (120021) 
andlRus zkowsk i fc BegehraaTT^OOoT) . In a similar approach, 
iMeverl l|200^ ~Droposed that in strong magnetic fields under 
radiation pressure, discs fragment into individual magnet- 
ically aligned columns of cool disc gas and that radiation 
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escapes through the gaps between them at super-Eddingtori 
luminosity. 

So far most investigations assumed a global limit in 
analogy to the classical Eddington limit to be valid also 
in th e context of an accretion disc. However. iFukuel (2000, 
2004) questioned this assumption for the first time quan- 
titatively and calculated a local Eddington-type limit for 
radiation pressure dominated geometrically thin discs as 
well as for slim discs by applying self - simila r solutions, 
based on the work of IWatarai fc Fukud dl999h . He found 
that inside a certain critical radius s cr , the disc remains 
in an Eddington-critical state with decreasing local accre- 
tion rates M = Mi npu t ■ (s/s cr ), independent of the disc 
model. The critical radius scales with the mass input rate, 
s cr ~ 1.95(Afin pu t/MB) • s g , with Me being the classical Ed- 
dington rate and s g the gravitational radius of the central 
black hole. 

In the present paper, we calculate the analog to the Ed- 
dington limit (resp., the critical accretion rate) for a classi- 
cal thin disc model of the lShakura fc S unyaevl Jl973l) variety 
and for a slim disc model of the lAbTamo"wfe^t^rt] jl98al 
type. In both descriptions, a standard a- viscosity is applied, 
therefore restricting the models to the non-selfgravitating 
case. We relinquish the above approximations (1), (2), (3), 
(i), and (ii), but keep the approximation of non-relativistic 
stationarity. 

In Sect. |21 we provide the models and derive the equa- 
tions for the critical accretion rates for both disc types and 
directions (vertical, radial) separately. In Sect. [3] the results 
are presented and compared to the classical Eddington ac- 
cretion rate. Based upon that, a locally Eddington limited 
disc model is proposed. Sectionals dedicated to the discus- 
sion of our results, and Sect.|5]to an outlook on the possible 
effects on the evolution of accretion discs. 



2 MODELS 

The original physical reasoning for the Eddington limit is an 
equilibrium of the gravitational force F s and the radiative 
force _F r . Contrary to a spherically symmetric body like a 
star, an accretion disc has to be treated as at least a two- 
dimensional object. Therefore, one can define two different 
"Eddington limits", corresponding to the equilibria in the 
vertical (z) and radial (equatorial, s) direction: 



F&{s) = Fy>(s)-F^(s) + 



F}*\s)-F^(s) 







= 



For a rotating viscous disc, further contributions to the total 
force have to be included. Moreover, the disc Eddington limit 
is no longer independent of the distance from the center of 
the object, as is true for the stellar case, leading to local, 
position-dependent Eddington limits. 

This section is divided into three parts: The first part 
introduces the general set-up which applies to thin and slim 
discs likewise, while the second and third part focus on thin 
and slim disc models, respectively. 

1 The classical Eddington accretion rate Mg is derived by equat- 
ing the disc luminosity with the Eddington limit and deriving a 
mass flow rate from it (see Sect. 12.11 



2.1 General set-up 

In this contribution, we use a cylindrical coordinate system 
{s, <p, z} with the distance r to the origin, r 2 = s 2 + z 2 . 

The critical accretion rate at a given position s is de- 
noted by Merit = M C rit (s), while the classical Eddington 
accretion rate is identified as Mb. 

If one naively identifies the classical Eddington lumi- 
nosity 



4-7TC 



GM 



(1) 



(k os = 0.4cm 2 g 1 for Thomson-scattering) with the total 
energy output of an accretion disc of 



GMM e 



(2) 



Me 



(3) 



one may derive a quasi-classical Eddington accretion rate 

4 si 

— • 7TC . 

Here, Sj denotes the disc's inner radius; the factor rj takes 
care of the coupling between the innermost disc regions and 
the accreting object. In this paper, we use r\ = 1/2. As long 
as 7] is of order unity, our conclusions do not depend on its 
exact value. 

We assume an optically thick situation for all our calcu- 
lations. This, together with the assumption that selfgravity 
plays no role, has to be justified afterwards. 

In contrast to the classical approach, we allow for opac- 
ity sources other than Thomson scattering and use an ana- 
lytic interpolation formula (Gail, pr iv.comm.; for a similar 
approach, see also lBell fc Lir] J1994D ): 



1 

^ in 



(3 000K) 1 



(3 000K) 1 
1 



■ T 1 



+ K mol 



+ 4- 



+ 



dust, evap 

The individual contributors axe approximated by 

Ki = «o,i ■ T KT - { ■ p Kp ' i 



1/4 
1/4 

(4) 
(5) 



and are compiled in TableQ 2 However, K- m is only meant as 
a handy interpolation. In particular in the transition regions 
between individual contributions, one has to be aware of this 
and should refrain from a physical over-interpretation there. 

We also take into account gas pressure, giving the total 
pressure as 

2pfc B T 4aT 4 

Ptot = Pgas + Prad = \ 7, • (6) 

mH 3c 

A parameter describes the contribution of gas to the total 
pressure in a standard way: 



Pgas 
Ptot 



1 - P 2aT 3 m H 



3cpfcE 



(7) 



Slight differen ces, for in stance for Ki CCj ovap by a factor of 2, 
between iBell fc Linl J1994 . Table 3) and our values, are due to 
somewhat different ways of combining the individual contributors. 
For the purpose of the present investigation, these differences are 
of minor importance 



© * RAS, MNRAS 000, HE] 



On the Eddington limit in accretion discs 3 

Table 1. Interpolation of the opacity: set of parameters (in cgs-units) 



Contributor i Symbol 



Dust with ice mantles 


^icc 


2.0 ■ 10~ 4 





2 


Evaporation of ice 


^ice, evap 


1.0 ■ 10 16 





-7 


Dust particles 


^dust 


1.0 ■ lO" 1 





1/2 


Evaporation of dust particles 


^dust, evap 


2.0 ■ 10 S1 


1 


-24 


Molecules 


^mol 


1.0 ■ 10~ 8 


2/3 


3 


Negative hydrogen ion 


K H~ 


1.0 ■ 10~ 36 


1/3 


10 


Bound-free and free-free-transitions 


^atom 


1.5 ■ 10 20 


1 


-5/2 


Electron scattering 


K - 
e 


0.348 









A main point of interest is the total luminosity of an Edding- 
ton limited accretion disc. We calculate the disc's luminosity 
by 

L disc = 2-27r- f ° dsaT* s , (8) 

J Si 

where T c s stands for the effective temperature (surface tem- 
perature) of the accretion disc at radius s, which is related 
to the central temperature T c by the optical depth r = Kph, 

°Tcs = ^crT c 4 , (9) 
or 

replacing the vertical stratification by a one-zone approxi- 
mation. 

Since our models take advantage of the a-description 
for the viscosity, the mass of the disc has to be 
small compared to the mass M of the central object 
iDuschl. Strittmatter fc Biermann 2000). For later compar- 
ison with the central mass, we introduce the disc's mass 

A/disc = / Z2n sds (10) 

with the surface density E = 2ph. The height h of the disc is 
defined by the vertical hydrostatic equilibrium, p/h = pg z , 
where the mass density p and the pressure p are measured 
in the disc plane (z = 0) and the gravitational acceleration 
g z is taken from the disc surface (z = h). 



while it is accelerated outwards by the radiative force 

F M = !™L aT * Ii (12) 
c 

according to the Stefan-Boltzmann law. Using 

4 3 GMM 
aT «* = ^^^ f (13) 

( z} ( z) 

and equating Fg and iy > we obtain an equation for the 
local critical accretion rate: 

/ \ 87T c h .„ 
Merits = — -•- 14 
3k/ 

Here, / is a function describing the radially inner boundary 
condition. A torque-free condition, for instance, leads to / = 

1 - yj Si/S. 

Equation 11411 changes the quality of the radial mass 
flow rate. While in standard disc theory, it is a parameter, 
it now becomes a local solution of our criticality assumption. 

By means of the remaining disc equations and J7J, we 
solve this set of equations for j3 and finally end up with 
P — 3/4, independent of a specific choice of the opacity n. 
This leads to 

M clit {s) = — —— § =75., (15) 

/ 3 3Va m J/ 2 Q l/8 K 9/8 (T l/8 n 7/8 

where S1k stands for the Keplerian angular velocity. Con- 
trary to the solution for /3, the critical accretion rate depends 
on the opacity. 



2.2 Geometrically thin accretion discs 



2.2.2 Radial direction 



Geometrically thin models have been the most popular de- 
scription of accretion discs in the last 30 years. Based upon 
the assumption of a negligible height h of the disc relative 
to the radius r and neglecting advective cooling, this model 
is simple but powerful for many applications. We derive the 
critical accretion rate separately for the vertical and radial 
direction by assuming a balance of force between the out- 
dragging radiative force F T and the attracting forces, like 
the gravitational force F s . 



2.2.1 Vertical direction 

In this case, we consider a particle with mass win at the 
surface h of the disc. This particle is attracted by the vertical 
component of the gravitational force 

GMniH h 



The calculation of the critical accretion rate given by a bal- 
ance of radial forces is less straight-forward then for the 
vertical direction. We consider a particle in the equatorial 
plane of the disc, at z = 0. So, the attracting gravitational 
force is given by 

{s) _GMrm_ {w) 



(11) 



In contrast to our previous derivation, the particle is located 
inside the disc and no longer at its surface. Assuming our 
particle being located at a specific position s in an opti- 
cally thick medium, additionally to the out-dragging radia- 
tion originating from the inner shell at s — ds we have to take 
into account also the radiation of the outer shell at s + ds 
pointing inwards. In brief, we have to apply the radial dif- 
ference of the radiative force instead of its absolute value. 
Furthermore, we have to consider the centrifugal force due 
to azimuthal motion and, coupled with it, the viscous forces. 
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Putting all together, this leads us to the radial component 
of the Navier-Stokes equation, which is given by 



2.3.1 Vertical direction 



P 



dpto 



+rriHV- 



d 2 v a 



Vs 

s 1 



TTlH dv 2 s 

GMmn 



= 0. 



(17) 



By means of the disc equations, the condition of a vanishing 
total force yields a differential equation for /3. In the Kep- 
lerian case, the centrifugal and gravitational force (i.e., the 
first and the last term of <\17l ) cancel each other, leaving the 
remainder of HI 711 to be solved. 

Independent of the chosen boundary values, its solution 
leads to supersonic radial motion, \v s \ > c s = (ptot/p) 1 ^ 2 , 
which contradicts the thin disc model. Hence, we investi- 
gate the corresponding value of the accretion rate under the 
condition v s = Ac s , \A\ < 1, which can subsequently be 
expressed as 



M cr it (s) = 



32ttc A 
3k a 



(1 - 0) 



(18) 



2.3 Slim accretion discs 



lAbramowicz et all jl98cf) extended the model of supercrit- 
ical accretion discs by I Jaroszvhski et alJ il98Cf) to slim ac- 
cretion discs with h ;$ s. Since the height h is no longer 
negligible compared to the equatorial distance s, the radial 
momentum equation has to be solved in a more precise way, 
leading to a non-Keplerian rotation of the disc. Furthermore, 
beside of the standard cooling by radiation, also cooling by 
advective mass flows becomes important, requiring a more 
detailed treatment of the energy (transport) equation. 

As the vertical disc structure is known only approxi- 
mately, the authors used three correcting factors (Bi, B 2 
and B3), which are all of the order unity, to correct the 
results of the vertical integration: 3 

• B\ takes into account that the height of the disc is 
no longer negligible compared to its radius, (s 2 + /i 2 ) 3//2 ~ 
Bi s 3 . In the model, we use a Newtonian potential; combined 
with the assumption of hydrostatic equilibrium, this leads 
to 



Ptot 

p 



GM ■ h 

(s 2 + h 2 ) 3 / 2 



GM ■ h 
Si -s 3 



Bi 



(19) 



• B2 is a measure of how efficiently the viscously dissi- 
pated energy is converted into radiation (rather than radi- 
ally advected), 



P 



m H 



(20) 



• B3 measures the efficiency of the advective energy 
transport in the energy transport equation and as such is 
part of modeling the advective processes. 



3 T hese parameters where int roduced originally 
by Paczvriski & Bisnovatvi-Koean (1981). However, the no- 
tation of the parameters differs slight ly between these tw o 
publications. In this paper, we follow lAbramowicz et all ll9. 



disc equations is given in 
p. 648). Starting from there 



The resulting s et of 
lAbramowicz et, ail l|l988 
and using the parameter (3 from Q, we derive two equa- 
tions for the three unknowns /3, h and I, where / is the local 
specific angular momentum. The missing third equation is 
again given by the Eddington condition of a vanishing total 
force: 



A*) 



= f! z) + r (z) = 



With 



GMrriH 



( s z + h 2 ) Vs 2 + h 2 
in terms of the slim disc description, and 



GMmuh 
B lS :i 



(21) 



(22) 



F, 



B-i frad/l = (1 

P 

B 2 GMmuh 



0)B 2 



hp 



(23) 



(24) 



the Eddington condition becomes 

Also in this model, the ratio of gas pressure to total pres- 
sure remains constant and is determined by the choice of 
the parameter B 2 . The critical accretion rate as well as the 
other disc quantities can be calculated by solving the two 
(coupled) differential equations numerically (Appendix 1X1. 

2.3.2 Radial direction 

We find that similar to the thin disc case presented in 
Sect. l2.2"!?l the radial Eddington limit in slim accretion discs 
leads to supersonic radial velocities, \v s \ > c s , which cannot 
be the case for (stationary) slim discs. Therefore, we follow 
the same idea as in Sect. |2~2~^1 and require v s = Ac B with 
L4I < 1. In terms of the slim disc model, we obtain 



The condition v s 



B x B 2 {l-l ) 
Ac s leads to 



h = B 2 V5r.4.(l-i<>',>-. 

a \Ik Ik, 



(25) 



(26) 



Here, Ik and lo are the local Keplerian specific angular mo- 
mentum and the specific angular momentum of the disc's 
innermost orbit, respectively. 

Starting from the radial momentum equation and the 
energy equation for slim discs, we can derive two additional 
differential equations for the unknowns h, f3 and I. Together 
with 12611 . the system is closed and can be solved numeri- 
cally. The set of equations is shown in appendix iBl 



3 RESULTS 

As can be seen from the above equations, the critical accre- 
tion rates depend on the distance s from the central object. 
We wish to emphasize that these accretion rates represent 
the local maximum value of mass which can be transferred 
inwards at a specific distance from the center in a stationary 
accretion state. 
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thin, vertical, el. scat, 
thin, vertical, int. opac. 
thin, radial, el. scat./int. opac. 



g/Sj 

Figure 1. Relative height h/s for the vertical Eddington limit 
and for the radial limit in standard o-discs 



M crit /M B 
o.oi 



thin, vertical, el. scat, 
thin, vertical, int. opac. 
thin, radial, el. scat., A=l/100 
thin, radial, int. opac, A=l/100 



1000 



10000 



s/sj 



Figure 3. Critical accretion rate M cr ; t in units of the classical 
Eddington rate Me for the vertical Eddington limit and for the 
radial limit in standard a-discs 



thin, vertical, el. scat./int. opac. 

thin, radial, el. scat, 
thin, radial, int. opac. 



s/sj 

Figure 2. /3 = p gas /ptot for the vertical Eddington limit and for 
the radial limit in standard a-discs 



For the following discussion of our results, we define a 
reference model: The central object is assumed to be a stellar 
black hole with M = IOMq, resulting in an inner radius of 
si = QGM/c 2 = 8.86 • 10 6 cm. The outer radius s is set to 
5 ■ 10 4 Si. In this environment, the classical Eddington rate 
is M E = l,67-10 19 gs _1 = 2.65 • 10~ 7 M© a -1 . The viscosity 
parameter a is set to 0.1. A torque-free boundary condition 
is chosen at the inner disc radius s = si: f = 1 — \J Si/s. 

In the following figures, the curves are labeled with 
the corresponding disc model (thin, slim), local Eddington 
limit (vertical, radial) and opacity (el. scat.: k — k cs , 
int . opac . : k = Km). 



3.1 Thin accretion discs 

FiguresQ£|]present the results for the vertical and the radial 
direction - in the latter case, the parameter A is set to 0.01, 
corresponding to h/s ~ 0.1 as an upper limit for thin discs. 
This choice is somewhat arbitrary, but provides a resonable 
upper limit for the radial velocities in thin discs. According 
to 1181 1. the allowed mass accretion rate scales linearly with 
the parameter A. 

Figure Q shows that the height of Eddington limited 



discs decreases steadily with the radius s, while it remains 
constant for a large range in the radial limit. This is a direct 
consequence of the gas-to-total-pressure ratio /3, which is 
constant only for the vertical Eddington limit (Fig.[5Jl. In the 
radial limit, /3 decreases with s by two orders of magnitude, 
causing the radiation pressure to exceed the gas pressure 
and to expand the disc vertically. This leads to higher sur- 
face densities and radial velocities, and therefore to higher 
critical accretion rates, as it can be seen from Fig. [3] Further- 
more, we observe a notable influence of the opacity on the 
results: For k os , the critical accretion rates are much higher 
than for Ki n , reaching a value of approximately the classical 
Eddington rate Me in the inner part of the disc. According 
to 112H and the fact that in thin discs all energy released by 
viscous processes is transported by radiation, the radiation 
force depends linearly on k, and Af cr it oc K , see 1141 . 

For some parameter ranges we find rather large differ- 
ences between the Eddington limits derived following the 
standard approach of taking into account electron scatter- 
ing as the only absorbing process («es), and those which al- 
low for additional absorption processes (Km). In the present 
calculations, bound-free and free-free processes (contributor 
K a tom in 0) dominate the opacity jq n and alter it signifi- 
cantly. The resulting differences in the Eddington limits are 
of relevance and emphasize the importance of taking into ac- 
count absorption processes other than Thomson scattering 
in the disc case, i.e., of using consistent opacities. 

The changes in the inner few s; is due to the boundary 
condition / at the inner radius: For s — * Si the boundary 
function / tends to zero. As M cr it scales with / _1 (see 11511 ), 
the critical accretion rate formally diverges. This, however, 
has to be taken as a technical artefact: In reality, the stress at 
the inner boundary will always be nonzero, leading to finite 
values for Merit. However, / may be still be small and there- 
fore theoretically allowing for higher accretion rates than at 
the minimum at s m i n ~ 2si. As long as no vertical mass in- 
flows to the disc are present, the real accretion rate M will 
stay at its minimum value M cr it(s m in) and the disc will be in 
a sub-Eddington state in terms of our local Eddington-limit 

for Si < s < S m i n . 

For the vertical Eddington limit and also for the ra- 
dial limit, the condition of a non-selfgravitating disc holds, 
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slim, vertical, el. scat./int. opac, h/s<=2.00 
slim, vertical, el. scat./int. opac, h/s<=1.00 
slim, radial, el. scat./int. opac. 



h/s 



slim, vertical, el. scat./int. opac, h/s<=2.00 
slim, vertical, el. scat./int. opac, h/s<=1.00 
slim, radial, el. scat./int. opac. 



s/Sj 

Figure 4. Angular momentum I in units of the Keplerian angular 
momentum Zj< f° r the vertical Eddington limit and for the radial 
limit in slim discs 



s/sj 

Figure 5. Relative height h/s for the vertical Eddington limit 
and for the radial limit in slim discs 



Mdisc -C M. Also the assumption of optical thickness is ful- 
filled for both K es and Ki n , 



3.2 Slim accretion discs 

For the model calculations, we choose the three parameters 
Bi, B2, and B3 as follows: 

• The slim disc model is valid for h J$ s, i. e. (s 2 + /i 2 ) 3 / 2 = 
Bis 3 J$ 2 3//2 s 3 , making B\ — 2 a reasonable choice. 

• For the vertical Eddington limit, j3 = 1 — I/-B2 • In order 
to compare our results directly to thin accretion discs with 
(3 = 3/4, we set B 2 = 4. 

• We follow iPaczvriski fc Bisnovatvi-Koganl jl98lT) and 
set B 3 = 1/2. 

In order to solve the sets of differential equations, we 
have to define boundary conditions. Due to numerical rea- 
sons, this has to be done at the outer boundary. We set 
l(s ) = Ik(s ) and /3(s ) = 3/4 for the radial limit. With 
A — 1 /30 (for slim discs, radial velocities are usually higher 
than for thin discs) and 1261 . this leads to h(s Q ) = 2s . Also 
here, the choice of A is somewhat arbitrary. Since the depen- 
dency of the allowed mass accretion rates on A cannot be 
seen directly from the equations, we also compute the case 
A = 1/20. For the vertical limit, we set l(s a ) = ^k(so) and 
h(s ) — 2s a . For comparison, we also compute the vertical 
limit with h(s ) — s . The results are shown in Figs.^HZI 

As it can be seen from Fig. |3] the angular momentum 
drops with respect to the local Keplerian value Zk(s) for 
both the radial and the vertical limit. For the latter one, the 
decrease is extremely steep, resulting in l(s) — (l + e)'o, £ <S 
1. For s — * si, the angular momentum reaches its Keplerian 
value in the vertical limit, while it overshoots this limit at 
s ~ si in the radial limit. 

Figure |S] shows the behaviour of the disc's relative 
height h/s for local equilibria of forces in the vertical and 
radial direction, respectively. While this ratio scales with s 
in the vertical limit, it adopts a constant value about unity 



slim, vertical, el. scat/hit. opac, h/s<=2.00/l .00 
slim, radial, el. scat./int. opac. 



s/si 

Figure 6. /3 = p gas /ptot for the vertical Eddington limit and for 
the radial limit in slim discs 



M crit /M E 

le+008 



slim, vertical, el. scat./int. opac, h/s<=2.00 
slim, vertical, el. scat., h/s<=l .00 
slim, vertical, int. opac, h/s<=1.00 
slim, radial, el. scat./int. opac, A=l/30 
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Figure 7. Critical accretion rate M cr ;t in units of the classical 
Eddington rate A/e for the vertical Eddington limit and for the 
radial limit in slim discs 
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Figure 8. re es and rei n for slim discs in the radial limit 



in the radial limit. In analogy to thin discs, this corresponds 
to a decrease of f3 in the latter case for smaller radii (Fig. |SJ 
and to higher accretion rates with respect to the vertical 
limit (Fig. 0. The critical accretion rates exceed the limits 
for thin discs (cf. Sect. I3.l"l significantly. 

We want to point out that the results are almost in- 
dependent of the opacity, although the absolute values of k 
differ considerably (Fig. [HJ. This becomes plausible wh en re- 
garding the energy equation jAbramowicz et al Jl988l . (12)): 
For high accretion rates, the radiative processes are negligi- 
ble compared to the advective flows. 

In all cases, the assumption of large optical thickness is 
fulfilled. By computing realistic disc models (see Sect. l3.3l for 
a more detailed discussion), we find that also the condition 
of a negligible disc mass is fulfilled as long as the mass input 
rate at the outer boundary M(s ) stays at a reasonable value 
below 10 5 Af E . 



3.3 An Eddington limited thin disc 

We have seen in Sect. |3~T1 that the Eddington limit poses 
strong limitations on the mass inflow rate in classical thin 
accretion discs. Contrary to the above sections where we 
displayed the disc properties for critical accretion rates at 
every radius, we now construct a stationary thin disc model 
by taking into account the local Eddington limit only where 
needed. 

Let us assume an external mass input rate M at the 
outer boundary s . As the critical accretions rates in the 
outer part of a thin a-disc become quite large, we expect 
the disc there to be sub-critical for any reasonable M . 
But, at a certain radius s cr i t - which can be calculated di- 
rectly from 1151 - the external mass input rate may become 
equal to the local Eddington limit Merit (scrit). As M CT it{s) 
is monotonically decreasing (for the influence of the inner 
boundary condition, see the discussion above), the disc will 
remain in a critical state for all s ^ s cr i t when assuming that 
mass outflows guarantee the condition M(s) — M cr it(s) for 

S ^ Scrit • 

This leads to discrepancies in the disc model, as the a- 
disc equations have been derived under the condition M = 
const. This requires the application of non -conservative disc 
models as presented bv iLipunoval lll999l) . As can be seen 



M/M] 




Figure 9. Accretion rate for a thin disc, M = 10M^. Solid/ 
dashed line: with/without consideration for the local (vertical) 
Eddington limit 



h/s 




Figure 10. Relative disc height h/s for a thin disc with M a = 
lOMg. Solid/ dashed line: with/without consideration for the lo- 
cal (vertical) Eddington limit 



from her results, the disc equations differ only slightly from 
that of a conservative disc model with constant accretion 
rate. This is also what we expect from comparing the time 
scales: The disc will adjust to changes in the accretion rate in 
about the hydrostatic timescale fhyd = h/c s , which itself is 
much shorter than the viscous timescale i v j s = s /v s , as long 
as the disc remains thin. Therefore, we can approximate an 
Eddington limited thin disc with a quasi-static model. 

Fi eures 1 1 1 show the results for such a thin disc with 
M = IOMe. For simplicity, we adopted simple electron 
scattering throughout the whole disc and ignored the in- 
ner 2si because of the influence of the boundary condition. 
In this case, s cr it ~ 750si. As can be seen from these figures, 
the mass input rate to the central object drops to about 
10 _2 Me in the range [2si;s cr it]. Also, the ratio h/s follows 
this trend, in contrast to the pseudo thin disc solution shown 
in dashed lines with M(s) = M a = const. The disc inflation 
in the latter conservative case corresponds to a transition 
into a radiation pressure supported disc and violates the as- 
sumptions of geometrical thinness. On the other hand, this 
condition of negligible relative height of the disc holds for 
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Figure 11. j3 = p ga s/ptot for a thin disc, M a = IOMe- Solid/ 
dashed line: with/without consideration for the local (vertical) 
Eddington limit 

an Eddington limited thin disc, which is mainly supported 
by gas pressure. 

The total disc luminosity, given by ©, becomes 0.15Le 
for the Eddington limited case and IOLe for the pseudo thin 
disc case. The total disc masses 1110ft almost equal each other, 
A^disc ~ 10 _6 M (M: mass of the central object), reflecting 
the fact of a simple "puffing up" of the pseudo thin disc. 



4 DISCUSSION 

In order to interpret our results correctly, we have to keep in 
mind that the disc quantities presented in Sect . 13 . ll and 13 . 2l 
correspond to discs at their respective local vertical and ra- 
dial Eddington limits. The critical accretion rates in the 
outer regions are enormous. However, such high mass in- 
put rates are only formal solutions as they will never be 
reached in real systems. The allowed mass accretion rates 
are decreasing constantly with smaller radii, confirming our 
expectations that the Eddington limit becomes important 
only in the inner part of the disc. 

In all our models, a strong correlation between the rel- 
ative height ft/s of the disc and the critical accretion rate 
Merit is present. For a constant ratio of gas pressure to total 
pressure /3, the relation h/s scales like s 0,6 . From the anal- 
ysis of the radial limit, we know that with decreasing /3 for 
a constant radius, the disc is blown up in vertical direction, 
i.e. h/s increases. Associated with this development, but 
less effective, Merit increases too. We have also seen that for 
sufficiently high M and h/s, the angular momentum in the 
disc differs significantly from its Keplerian value over large 
parts of the slim disc, allowing advective flows to dominate 
the energy transport and to rise the critical accretion rates. 

We can now draw a more realistic scenario for an ac- 
cretion disc system: Let us assume a given reasonable mass 
input rate M at the outer boundary of the disc. Doing so, 
we can certainly assume the disc to be thin in its outer parts. 
For smaller radii, the ratio h/s will begin to increase as long 
as the vertical Eddington limit is not violated and therefore 
no outflows occur, see Fig. 1101 Further inwards, two compet- 
ing effects now come into play: On the one hand, outflows 
may occur and help the disc to stay in a subcritical state; 



on the other hand, the thin disc may pass into a slim disc 
with deviations from Keplerian rotation, advective energy 
transport and a changing ratio f3. In reality, both effects 
will be in concurrence. To answer this question, extended 
time- dependent disc models have to be investigated. 



5 CONCLUSION 

Our work shows that the quasi-spherical treatment of accre- 
tion discs by applying the classical Eddington limit does not 
hold as a good approximation. The allowed mass accretion 
rates differ significantly from the stellar case and also de- 
pend strongly on the disc model, while the disc luminosity 
only weakly exceeds the classical Eddington value. 

We have seen that the inner domain of an accretion 
disc holds the key position for the evolution and also for the 
observational appearance of highly accreting systems: The 
final growth rate of the central object, the total amount 
of matter that has to be expelled from the disc and also 
most its luminosity are determine d by it. Numerous ra dia- 
tion hydrodynamic simulations (see O hsuga et alj 1120051) for 
one example) show that outflows in the middle part of the 
disc get stuck and may therefore be driven back to the disc, 
while mass ejections in the inner part lead to high velocities 
and gravitationally unbound streams of previously accreted 
material. 

It is broadly accepted that photon tr apping effects, 
which are not included in the slim disc model llOhsuga et alJ 
[2002j), alter the emerging luminosity and therefore may 
change our results in a way that even higher critical ac- 
cretion rates are possible. Also, the effect of further energy 
transport mechanisms, like convection and heat conduction, 
and non-stationarity have to be investigated in a more so- 
phisticated discussion. 
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APPENDIX A: VERTICAL EDDINGTON LIMIT 
IN SLIM ACCRETION DISCS 

The parameterization of the pressure allows us to reduce all 
thermodynamics to the parameter /3: 



Ptot 
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We use the resu lt <|24^. 8 = cons t ., tog ether with equa- 
tions (3)-(13) of lAbramowicz et alJ l)l988T ). More precisely, 
we substitute all unk nowns except of h and I in the radial 
momentum equa tion jAbramowicz et al.|ll988t (8)) and the 
energy equation ijAbramowic^t^L^^sM (12)) by using the 
remaining disc equations and the Eddington condition 1248 . 
Note that (7) contains an error in the second term in paren- 
theses: The factor (4 — 8) has to be replaced by (4 — 3/3) 
llJaroszvriski et al"lll980l) . 

This leads us to two differential equations for the disc 
height h and the angular momentum I, namely 



Herein, Ik (resp. S7k) stands for Keplerian azimuthal motion 
and lo is given by the boundary condition of a vanishing 
viscous torque at the inner boundary si of the accretion disc, 
= 'k( s 0- The critical accretion rate can be calculated 
from the solutions of IA4I and IA5t by 
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APPENDIX B: RADIAL EDDINGTON LIMIT 
IN SLIM ACCRETION DISCS 

Again, we use the /3-description for the disc's thermodynam- 
ics, which yields to 1A1I - IA3I . The missing equation to be 
able to determine M cr i t from the system is now given by 
v s = Ac s instead of the Eddington condition 12 H . From the 
disc equations and from the first law of thermodynamics, we 
already derived 12611 : 



03 



As we have seen in Sect. 12.2 - !?! we cannot assume = const, 
in this case. The radial momentum equation and the energy 
equation lead to two differential equations for h, I and 8, 
namely 
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The critical accretion rate is again given by 1A6|I . Together 
with (I26H . this system can be solved numerically. 
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